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Abstract — We analyze pattern formation on a network of 
cells where each cell inhibits its neighbors through cell-to-cell 
contact signaling. The network is modeled as an interconnection 
of identical dynamical subsystems each of which represents the 
signaling reactions in a cell. We search for steady state patterns 
by partitioning the graph vertices into disjoint classes, where 
the cells in the same class have the same final fate. To prove 
the existence of steady states with this structure, we use results 
from monotone systems theory. Finally, we analyze the stability 
of these patterns with a block decomposition based on the graph 
partition. 

I. Introduction 

Spatial patterning plays a crucial role in multicellular 
developmental processes [1], [2]. The majority of theoretical 
results on pattern formation rely on diffusion-driven 
instabilities, proposed by Turing [3] and further studied by 
other authors [4], [5], [6]. Although some early events in 
developmental biology employ diffusible signals, most of the 
patterning that leads to segmentation and fate-specification 
relies on contact-mediated signals. Lateral inhibition [7] 
is a cell-to-cell contact signaling mechanism that leads 
neighboring cells to compete and diverge to contrasting 
states of differentiation. An example of lateral inhibition 
is the Notch pathway, where neighboring cells compete 
through the interaction between Notch receptors and Delta 
hgands [8]. 

Dynamical models for the Notch mechanism have been 
analyzed in [8] and [9]. However, these either include 
analytical results for only two cells, or perform numerical 
simulations for larger networks. A broader dynamical model 
for lateral inhibition is proposed in [10], and results that 
are independent of the size of the network are presented. 
In this reference, the large-scale network is viewed as 
an interconnection of individual cells, each defined by an 
input-output model. The contact signaling is represented by 
an undirected graph, where each vertex is a cell, and a link 
between two vertices represents the contact between two 
cells. Results for the instability of the homogeneous steady 
state and the convergence of two level patterns for bipartite 
contact graphs are presented in [10]. 

In this paper, we use the model introduced in [10] and 
derive results for pattern formation on a general contact 
graph, recovering the results of [10] for bipartite graphs as a 
special case. Our main idea is to partition the graph vertices 



into disjoint classes, where the cells in the same class have 
the same final fate. We use algebraic properties of the graph 
and tools from monotone systems theory [11] to prove the 
existence of steady states that are patterned according to 
these partitions. Finally, we address the stability of these 
patterns by decomposing the system into two subsystems. 
The first describes the dynamics on an invariant subspace 
defined according to the partition; and the second describes 
the dynamics transversal to this subspace. 

A key property that each partition must satisfy is that 
the sum of the weights from one vertex in one class to 
those in another class is independent of the choice of 
vertex. Partitions with this property are called equitable 
and allow us to study a reduced model where all vertices 
in the same class have the same state. As examples of 
equitable partitions, we study bipartite graphs and graphs 
with symmetries. For symmetric graphs, we show that 
subgroups of the automorphism group of a graph can be 
used to identify equitable partitions. 

The idea of grouping the vertices of a network into 
classes of synchronized states has been explored in [12], 
[13]; however, no investigation about steady-states and 
their stability is pursued, and no biological application 
is addressed. Symmetry properties have been exploited 
in the dimension reduction and block decomposition of 
semidefinite program optimization problems, such as 
fastest mixing Markov chains on the graph [14], [15], and 
sum-of-squares [16]. Symmetry has also been related to 
controllability on controlled agreement problems [17]. 

In Section II, we define the model and introduce necessary 
graph theoretical concepts. In Section III, we present the 
main result of the paper, which provides conditions for the 
existence of steady-state patterns. In Sections IV, we apply 
the main results to graphs with symmetries, and present 
examples. Finally, in Section V, we present a decomposi- 
tion that is helpful for the stability analysis of steady-state 
patterns, and a small gain stability type criterion. 

II. Lateral Inhibition Model 

We represent the cell network by an undirected and 
connected graph Q — Q{V,E), where the set of vertices V 
represents a group of cells, and each edge e £ E represents a 



contact between two cells. The strength of the contact signal 
between cells i and j is defined by the nonnegative constant 
Wij. We let Wij — when i and j are not in contact 
and allow uneven weights to represent distinct contact signal 
strengths. This contact graph is undirected, i.e., W = {wij} 
is symmetric. 

Let N be the number of cells and define the scaled adjacency 
matrix P G R^^^ of ^ as: 

where the scaling factor is the node degree di = Tlj'^i.j- 
The definition of P implies that the matrix is nonnegative 
and row-stochastic, i.e., PIn — ^n, where l^v G 
denotes the vector of ones. The structure of P is identical to 
the transposed probability transition matrix of a reversible 
Markov Chain. Therefore, P has real valued eigenvalues 
and eigenvectors. 

Consider a network of identical cells i = 1, . . . , N whose 
dynamical model is given by: 

j — f[xi,Ui) 

\ Vi ^ h{xi) 

where Xi £ X d M" describes the state in cell i, G C M 
is an aggregate input from neighboring cells, and ?/i G 3^ C M 
represents the output of each cell that contributes to the input 
to adjacent cells. 

We represent the cell-to-cell interaction by 

u = Py (3) 

where P is the scaled adjacency matrix of the contact graph 
as in ([T]i, u := [ui...un]^, and y :— [yi . . .yN^'^ ■ This 
means that the input to each cell is a weighted average of 
the outputs of adjacent cells. 

We assume that /(•,•) and h{-) are continuously differen- 
tiable, and that for each constant input u* G U, system Q 
has a globally asymptotically stable steady-state 

X* := S{u*). (4) 

Furthermore, we assume that the map S :U ^ X and map 
T -.U ^y, defined as: 

T{-):=h{S{-)), (5) 

are continuously differentiable, and that T{-) is a positive, 
bounded, and decreasing function. The decreasing property 
of T is consistent with the lateral inhibition feature, since 
higher outputs in one cell lead to lower values in adjacent 
cells. 

Note that the steady states of the system ©-([Si) are given 
by Xi — S{ui) in which ui^...,un are solutions of the 
equation: 

u = PTjv(u), (6) 

where 

Tn{u) = [T{u{),...,T{un)V. (7) 



Since P is row-stochastic, (|6| admits the homogeneous 
solution Ui ~ u* for all i — 1,...,N, where u* is the 
unique fixed point of r(-), i.e, 

T{u) = u. (8) 

III. Identifying Steady State Patterns 

To identify nonhomogeneous steady states, we introduce 
the notion of equitable graph partitions. For a weighted and 
undirected graph Q{V.E) with scaled adjacency matrix P, 
a partition tt of the vertex set V into classes Oi , . . . , O,. is 
said to be equitable if there exist p^^, i,j = 1, . . . , r, such 
that 

^ Puv = Pij Vu G Oi. (9) 

This definition is a modification of [18, Section 9.3] which 
considers a partition based on the weights of the graph Wij 
instead of the scaled weights pij in ([U. 

The quotient matrix P G W"^"^ is formed by the entries 
Pjj. It is also a row-stochastic matrix, and its eigenvalues 
are a subset of the eigenvalues of P, as can be shown with 
a slight modification of [18, Thm9.3.3]. As we will further 
discuss, equitable partitions are easy to identify in bipartite 
graphs, and in graphs with symmetries. 

We search for nonhomogeneous solutions to (|6]l in which 
the entries corresponding to cells in the same class have the 
same value. This means that we examine the reduced set of 
equations 

z = PT,(z), (10) 

where P is the quotient matrix of the contact graph Q, and 
z . The patterns determined from the solutions of ( fTOl l 
are structured in such a way that all cells in the same class 
have the same fate, i.e, 

Ui — Zj for all i G Oj. (11) 

We now present a procedure to determine if ( fTOl l has a non- 
homogeneous solution. Define the reduced graph Q.^ to be a 
simple graph in which the vertex set is V = {Oi, . . . , 0^} 
and the edge set is 

E = {(0„ Oi) -.1^3. P^, 7^ or pj, ^ 0}. (12) 
Note that we omit self-loops in Q.^ even if 7^ 0. 

Assumption 3.1: The reduced graph Q-^ is bipartite. 

In the following theorem, we determine whether there 
exists a solution to the reduced set of equations ( fTOl i other 
than the homogeneous solution z* — u*lr. 

Theorem 3.2: Let tt be an equitable partition of the ver- 
tices of Q such that Assumption 13.11 holds. Let Vr be the 
eigenvector of P associated with the minimum eigenvalue 



Xr- If T{-) is positive, bounded, and decreasing, and if 
T'{u*) is such that 

\T'{u*)\Xr<-l, 
then there exists a solution of ( fTOl i other than 



Z* = U*lr. 

Proof: Consider the auxiHary dynamical system: 



(14) 
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:= F(z), z e K->o. 



Note that around the homogeneous steady state z* 
Jacobian matrix 

DF{z*) = -Ir + T'{u*)P 



from the fact that at z = 0, i = Pr(0) > 0, and at z = m, 
z < (since u > T{u)). The sets 

(13) Si = {u*l + R'^o)n[0,uY and 5*2 = (^t*l-i^?>o)n[0,^i]^ 

(19) 

are forward invariant. Therefore, we conclude from [11, 
Theorem 4.3.3], there exists an equilibrium point z ^ z* , 
and it satisfies ( fTol i. ■ 

Example 1: Checkerboard Patterns in Bipartite Graphs 
Suppose that the contact graph Q is bipartite, and choose Oi 
and O2 to be the partition such that every edge can only 
connect a vertex in Oi to a vertex in 02- Then, up to vertex 
relabeling, the scaled adjacency matrix of Q can be written 
as 

P12 
P21 



(15) 
, the 



(16) 



is a nonpositive matrix (since P is nonnegative, and from 
(O, T'{u*) < 0). 

We show that under a coordinate transformation the system 
is cooperative, see [11, Definition 3.1.3]. Following the 
bipartite property of Q^r in Assumption 13.11 we define a 
partition J C {1, . . . , r} and J' = {1, . . . , r}\J such that 
no two vertices in the same set are adjacent. Let ej = if 
j £ J and ej — 1 if j E J', and choose the transformation 
Rz to be 



i^ = d^ag{(-l)'^...,(-l)^^...,(-l)^'-}. 



(17) 



Since the reduced graph is bipartite, R^^PR — RPR is a 
matrix similar to P and all of its off-diagonal elements are 
nonpositive. In the new coordinates Rz, the Jacobian matrix 
in (fTST i becomes J :— R{DF{z*))R and has nonnegative 
off-diagonal elements. This means that the system is coop- 
erative. 

To prove the existence of a solution z ^ z* of (fTOb . we 
appeal to [11, Theorem 4.3.3] which stipulates that the 
largest real part of the eigenvalues of J (designated as s( J)) 
to be positive with associated eigenvector t; 3> {i.e., 
all elements are positive); and that there exists a bounded 
forward invariant set. 

First, note that J is a quasi-positive and irreducible ma- 
trix (this is because the reduced graph is connected, and 
T'{u*) ^ 0). Then, we know from [11, Corollary 4.3.2] that 
there exists an eigenvector v ^ such that Jv ~ s{J)v. For 
this case, the eigenvalues of J are all real and given by 



l + AfcT'(M*), fc = l,...,r, 



(18) 



where Xk are the eigenvalues of P. Therefore, s(J) ~ 
— 1 + T'{u*)Xr. From condition ( flj] ) we conclude that 
s(J) > with positive eigenvector v, and that Vr — Rv is 
an eigenvector of P associated with (i.e., PRv — XrRv). 
Second, since the transformed cooperative system is mono- 
tone with respect to the standard cone i?>o, we conclude 
that u*lr + i?>o and u*lr — i?>Q are forward invariant. 
Furthermore, since r( ) is bounded and decreasing (and 
T{u*) = u*), there exists an hypercube [0,u]'', with < 
u* < u which is also forward invariant. This can be seen 



P = 



(20) 



Since the rows of P12 (and also rows of P21) sum up to 
1, we conclude that tt, consisting of sets Oi and O2, is an 
equitable partition. Moreover, the reduced graph CJ^ is itself 
bipartite {i.e.. Assumption 13 . 1 1 holds), and matrix P is given 
by 

' 1 
1 



P 



(21) 



Since the eigenvalues/eigenvectors of P are Ai = 1 and 
A2 = — 1, the next Corollary follows. 

Corollary 3.3: Let Q be bipartite, and define a partition 
Oi C {1, . . . , A^} and O2 = {1, • ■ • , N}\Oi such that no 
two vertices in the same set are adjacent. Then, if 

\T'iu*)\ > 1, (22) 

there exists a steady state u — [ui . . .un] such that Ui = zi 
if i e Oi, and = Z2 if i G O2, and z\^ Z2^ u* . ■ 



The steady state defined by Corollary 13.31 results in a 
"checkerboard" pattern as in Figure|2|A), since it has distinct 
states for adjacent cells. 

IV. Graphs with Symmetries 

An important class of equitable partitions results from 
graph symmetries, which are formalized with the notion 
of graph automorphisms. For a weighted graph Q{V,E), 
an automorphism is a permutation g : V V such that 
if e E then also {gi,gj) G E and Wij = Wgi,gj, 

where gi denotes the image of vertex i under permutation 
g. The set of all automorphisms forms a group designated 
by automorphism group, Aut{G). A subset of a full 
automorphism group Aut{G) is called a subgroup if H is 
closed under composition and inverse. 
Let H he a subgroup of a full automorphism group Aut{Q). 
Then, the action of all permutations h E H forms a partition 
of the vertex set into orbits, Oi — {hi : h E H}, such 
that Oj = Oi for all j E Oi. Let r be the number of 
distinct orbits under the subgroup H, and relabel them as 
{Oi, . . . ,Or}. This orbit partition is equitable, because the 



E 



Pi'k 



is constant independently of the defined by an orbit partition. 



choice of i* e Oi. 

Since any subgroup of the full automorphism group of 
a graph leads to an equitable partition, we conclude by 
Theorem |32] that any orbit partition generated by a subgroup 
of Aut{Q) is a candidate for a pattern structured according 
to this partition. 

Example 2: Two-Dimensional Mesh 
Consider a two-dimensional mesh with wraparounds as 
in Figure [T] Since the graph is bipartite, an equi- 
table partition is given by the two disjoint subsets 
of vertices Oi = {1,3,6,8,9,11,14,16}, and O2 = 
{2, 4, 5, 7, 10, 12, 13, 15}. From Corollary [33] we know that 
a pattern with final value ui for all cells in Oi, and U2 for 
all cells in O2, with u\ ^ ^ u*, is a steady state of the 
network when \T'{u*)\ > 1; see Figure |2] (A). 



.10 !,n !,12 



Fig. 1. Graph representation for a two-dimensional mesli with wraparounds. 



The computation of automorphism groups, and the 
identification of the reduced order systems, becomes 
cumbersome as the size and symmetries of the graphs 
increase. However, these can be obtained from a computer 
algebra system with emphasis on computational group 
theory, such as GAP, [19]. 

Example 3: Two-dimensional Hexagonal Cyclic Lattice 
The number of distinct equitable partitions in a hexagonal 
lattice of cells is considerably large [20]. We use compu- 
tational algebra algorithms to find all the possible two-level 
equitable partitions obtained by automorphisms subgroups in 
a 6 X 6 cyclic lattice. Five distinct partitions, each with two 
classes, are plotted in Figure [3] 




(A) 



Fig. 2. Patterns for the graph in Figure (T). 

We next consider the automorphism subgroup that is gen- 
erated by a combination of two cell rotations in the horizontal 
direction, one cell rotation in the vertical direction, and one 
cell rotation in both vertical and horizontal directions. This 
subgroup leads to the orbits Oi = {1, 3, 5, 7, 10, 12, 14, 16}, 
O2 = {2,4,6,8,9,11,13,15}. The quotient matrix associ- 
ated with this partition is given by 



Pb = 



and has eigenvalues —1/2 and 1. Therefore, from Theorem 
13.21 a steady state state given by z = PT{z) exists if 
\T'{u*)\ > 2. 

In this example, the equitable partition obtained from the 
bipartite property of the contact graph can also be obtained 
by a subgroup of the automorphism group of the graph. 
However, this is in general not true; the four cell path is an 
example of a graph with a bipartite partition that cannot be 




Fig. 3. Five distinct partitions, each with two classes, in a 6 X 6 Two- 
Dimensional Cyclic Hexagonal Lattice obtained by symmetries on the 
contact graph. 



For these partitions, we have the following scaled adja- 
cency graph of the auxiliary system: 



1 

1 1 

2 2 



,Pc 



1 2 

! ! 

3 3 



and Pd = Pc, Pe = Pb- For each matrix we 
have the following smallest eigenvalue ~ —1/2, 

Xb = Xe = —1/3, and Ac ~ Xd =0. We thus conclude 
from Theorem 13.21 that the steady-state represented by 
pattern A in Figure |3] exists when |r'(w*)| > 2, and patterns 
B,E are steady-states when |T'(u*)| > 3. Theorem [32] is 
inconclusive for patterns C, D. 

Example 4: Soccerball Pattern on a Buckminsterfullerene 
Graph 

The next example addresses a larger graph, with 32 cells. It 
is motivated by the truncated icosahedron solid, also known 
as the Buckminsterfullerene [18], formed by 12 regular 
pentagonal faces, and 20 regular hexagonal faces, see Figure 
[4] In this case, we assume that each face is a vertex and that 




Determining stability for these steady states may become 
cumbersome for a large network of cells. To simplify 
this task, we decompose the system into an appropriate 
interconnection of lower order subsystems, and make the 
interconnection structure explicit. 




Fig. 4. Cell network and soccer ball pattern for the Buckminsterfullerene 
graph. 



two vertices are connected if the corresponding faces have a 
common edge. 

The full automorphism group leads to two orbits, one 
that consists of all the regular pentagon cells (Op), and the 
second orbit encloses all the regular hexagon cells (Oh)- The 
quotient matrix associated with the orbit partition is then 

1 

1/2 1/2 

This matrix has eigenvalues 1 and —1/2. Therefore, we 
conclude from Theorem 13.21 that a steady state as in Figure 
g] exists when > 2. 



P 



(23) 



Example 5: Nonbipartite nor Symmetric Equitable Parti- 
tion 

As discussed above, both bipartitions and automorphism sub- 
groups (symmetries) lead to equitable partitions. However, 
these are not the only cases that lead to equitable partitions. 
Consider the graph in Figure |5] with partition Ci ~ {3,6}, 
and 6*2 = {1,2,4,5,7,8}. This partition is equitable, but 
it does not result from an automorphism subgroup (for 
instance, there is no automorphism exchanging vertices 1 
and 4), and the graph is also not bipartite (due to the odd 
length cycles). 



1 



3 4 5 




P = 



1 

1/2 1/2 



Fig. 5. Example of an equitable partition that is neither bipartite nor a 
symmetry (Ci = {3, 6}, C2 = {1, 2, 4, 5, 7, 8}). 

The quotient matrix is the same as in Example 4, we 
conclude that a two level steady state pattern formed by 
cells Ci and C2 exists when |T'(\t*)| > 2. 



V. Stability Analysis of Patterns 

In the previous sections we discussed the existence of 
nonhomogeneous steady states for the cell network ©-(O. 



A. Network Decomposition 

Note that an equitable partition defines the subspace Xi = 
Xj for vertices i and j in the same class. This subspace is 
invariant for the full system (|2|i-(|3]l. Therefore, the steady 
states identified using an equitable partition of the contact 
graph lie on the corresponding invariant subspace. For a 
partition of dimension r, (Oi, . . . ,Or), the reduced order 
dynamics on this subspace consists of r subsystems as 
defined in (|2]i, coupled hy u = Py. 

Let the steady state be defined by Xi — S{ui), i = 1, . . . , N, 
where Ui — zj if i e Oj and [zi, . . . , Zr] is a solution of 
( [Tol l. The linearization at this steady state has the form 



x={A + BPC)x, 



(24) 



A 



(25) 



where A e t^Nti-kNu ^ block diagonal matrix where the 
?-th block is equal to Aj if i is in orbit Oj, and Aj £ R"^", 
is given by; 

dfix,u) 

Similarly, B e E^"x^ and C G M^""^" are block diagonal 
matrices as in A, with Bj = ^-^g^'"-* \ {xj.uj), and Cj = 

dh{x) I 
dx '^j' 

To decompose ( l24b into two subsystems, we select a repre- 
sentative vertex Vi for each class Oi. The set of r represen- 
tatives of each class defines the state of the subsystem on the 
invariant subspace. To see this, let Q be a matrix in R^^*", 
where qij = 1 if cell i is in class j, and qij — otherwise. 
Since the partition is equitable, we conclude that 



PQ - QP. 



(26) 



Letting T := [Q R] and choosing i? to be a matrix in 
^Nx{N-r) columns that, together with those of Q, form 
a basis for R^, we conclude that there exist matrices C and 
M such that 



P[Q R] = [Q R] 



P 




C 
M 



TP. 



(27) 



The matrix T is invertible and, thus, defines a similarity 
transformation from matrix P to P. Note that the upper 
left diagonal block of P is the matrix P, which describes 
the reduced order subsystem defined by the representative 
vertices. 

Next, we study a particular choice of the matrix R that 
gives a meaningful variable representation to the transverse 
subspace dynamics. Let the columns of R be given by 
standard vectors e^, defined as = Sij, j — 1,...,N; 
and further select the columns of R to be e^, i G 
{Oi\Vi, . . . ,Or\Vr}, in such a way that if i € Op\Vp, 
j e Oq\Vq, and p < q, then is in a column before ej. 



i.e., the column with non-zero entry i is placed before the 
column with entry j if vertex i is in a class with smaller 
index than the class of vertex j. 

For this choice of R we conclude from [21, Section 5.3] 
that the change of variables T (E) In, i-e., x — [T^^ ® In)x, 
leads to the decomposition of the linearized dynamics into 
the subsystem of the r representative cells Vi defining 
the invariant subspace, and the transverse subspace that is 
defined by the state difference between all the other cells in 
a class and their corresponding class representative: 



1, 



, ^ , where k G OAVj. 

(28) 

Therefore, the linearized system is decomposed into the 
representative subsystem Sr and the transverse subsystem 
Sd, 

Sr = {Ar,Br,Cr} and Sd ^ {Ad, Bd,Cd}. (29) 
To see this, note that 

i = {T-'^(Eln){A+BPC){T®La)i = {A+BPC)x, (30) 
where 

A T-'^ (g) I„AT (S In 

= diag{Ai,A2,...,Ar,Ai,... 

...,Ai,A2,...,A2,...,Ar,...,Ar} 

:= diag{AR,AD}, (31) 



with Ar e 



Ad e M(JV-r)«x(Ar-r)n^ matrices 



B and C are defined similarly to A. This means that to 
prove the stability of (|24] | we need to guarantee the stability 
of matrices 



Ar + BrPCr and Ad + BdMCd- 



(32) 



Stability certification of the steady state pattern is thus sim- 
plified due to the fact that P is typically much smaller than 
P, and M can be further decomposed by using a systematic 
approach based on the nested hierarchy of automorphism 
groups of the contact graph Q. 

B. A Small-Gain Criterion for Stability 

In this section, we provide a small-gain type condition for 
the stability of the steady state pattern around the solutions 
of (|6|l, based on the solution [zi , . . . , z^] of ( flOl l, which is 
defined by 



S{ui), i — 1, . . . ,N, where Ui = zj if i G Oj. 



(33) 



Consider again the linearization introduced in jZSl l. This 
describes a cell network with interconnection defined by 
u — Py and where each individual cell is given by the 
following linearized subsystem: 



A'x' 
C'x' 



B'u' 



(34) 



where — Aj, if i E Oj as in ( IZSl l, and similarly for 
and C\ see figure |6l Assume that each linearized subsystem 
Hi is observable and that A^ is Hurwitz. 



P * 



Fig. 6. Linearized System Interconnection. 

Since subsystems in the same class have identical models, 
they have identical £2-gains. Let 7,; denote the £2-gain of 
each subsystem in class i, and let F be a diagonal matrix 
with entries 



7i for j e 



(35) 



The following Proposition provides a small-gain criterion 
for the stability of the cell network around the steady state 
pattern. 

Theorem 5.1: Consider the network (|2]l-(l3]i. The steady 
state pattern defined by (|33] | is locally asymptotically stable 
if 

p{PT) < 1, (36) 

where F is as in ( [35] ). 

Proof: Since each linearized subsystem Hi, in ( [34] i, 
has bounded £2 gain, by the Bounded Real Lemma [22], we 
conclude that there exists a positive definite matrix Qi such 
that, for Vi{z) = z^QiZ, we have Vi{x'^,u'^) < jfu^'^u^ — 
y'^'^y'^. Let D = diag{di, . . . , dn}, where di is some positive 
constant, and let V{x) — J^i'^i^i- We then obtain. 



V{x) - }_^d,V^{x\u') 



(37) 



< U'^ DT^U - Y'^ DY (38) 
= Y'^{{TP)'^D{rP) - D)Y, (39) 

where U = [{u^f...{u^ff, Y = [{yY ■ ■ ■ iv^Vf , 
and the second equality follows from the fact that U = PY. 
If V{x) is negative semidefinite, we know from LaSalle's 
Invariance Principle [23], and for linear systems ([24]), that 
every trajectory converges to the unobservable subspace of 
the system. Furthermore, under the same assumption, and 
since each Hi is observable, we conclude that each trajectory 
must converge to the singleton {0}. Thus, the steady-state 
d33T l is locally asymptotically stable, if there exists a positive 
diagonal matrix D such that D — (TP)'^ D{rP) is positive 
definite. This is equivalent to the condition that / — TP be 
a M-matrix, see [25, Theorem 2]. To finalize the proof, note 
that the spectra of PT is the same as TP, and that the radial 



spectra assumption implies that / — TP is a nonsingular M- 
matrix, see [26, Definition 6.1.2]. ■ 



vertices in the same set are adjacent. As discussed before, in 
Example 1, the quotient matrix P is given by 



We now show that ( |36] | is equivalent to p(f F) < 1, where 

T = diag{'ji,...,jr}- (40) 

This result simplifies the verification of this small gain 
stability condition to the reduced system Sn in ( |29] l with 
interconnection matrix P. 

Lemma 5.2: Consider an equitable partition tt of Q. Then, 

p{PV) < 1 ^ p(PT) < 1, (41) 

where P is as in (|9]l, and F is as in ( l40b . 

Proof: To prove this statement we only need to show 

that 

p(PT) = piPT). (42) 

First, note that PT is a nonnegative irreducible matrix, by the 
Perron-Frobenius Theorem [26], we know that r = p{P^) > 
is an eigenvalue of PF with corresponding eigenvector 
iJ> 0. 

Claim: r > is also an eigenvalue of PT with corresponding 
eigenvector v such that entries Vi = Vj if i £ Oj. 
According to this claim, we know that u is a positive 
eigenvector Therefore, by citing again the Perron-Frobenius 
Theorem, and since PF is also a nonnegative irreducible ma- 
trix, we conclude that v has to be the eigenvector associated 
with eigenvalue r = p{PT). 

To prove the claim, note that matrix F is positive diagonal 
with repeated entries for vertices in the same class. There- 
fore, since the vertex partition is equitable for the scaled 
adjacency graph P, then it is also equitable when we consider 
a modified adjacency graph PF, i.e., 

^ Puv-fjj = Pijijj e Oi, (43) 

where p^j is as defined in Q. From this observation we see 
that this claim is a generalization of the Lifting Proposition in 
[14], which holds not only for partitions obtained through an 
automorphism subgroup but also for any equitable partition. 
The proof of the claim follows similarly to the proof of the 
Lifting Proposition, with matrices PF and PF as in i43[ . ■ 

Lemma 15.21 leads to a simplification of the stability con- 
dition ( |36] | in Theorem 15. II 

Corollary 5.3: Consider the network (|2])-(l3]l. The steady 
state pattern defined by (l33T l is locally asymptotically stable 
if 

p(pr) < 1. (44) 



C. Special Case: Bipartite Graph 

Let us consider the special case of a bipartite graph, with 
a partition tt consisting of two classes, chosen so that no two 



Therefore, the spectra of PF is ±\/7i72- The next result 
follows trivially from Corollary 15.31 

Corollary 5.4: Assume that the contact graph Q of the cell 
network (|2]i-(|3]l is bipartite, and that there exists a steady state 
u e such that Ui — zi if i G Oi and Wj = Z2 if i G O2, 
with 2i 7^ ^2 7^ z*, as in Corollary 13.31 Then, the steady 
state solution is locally asymptotically stable if 

7172 < 1, (46) 

where 71 and 72 are the £2 gains of the linearized subsys- 
tems around zi and Z2, respectively. ■ 

In the particular case where the £2 -gain is given by the 
dc-gain, 

7, = -C^AT'B' = -T'(z,), (47) 

we see that the local asymptotic stability condition in ( l46b 
reduces to 

r'(zi)T'(z2) < 1. (48) 

The £2 -gain is indeed equal to the dc-gain in the particular 
case where each subsystem (|2]l is input-output monotone 
[27], as assumed in [28]. We have thus recovered Theorem 
2 in [28] which used ( l48b to prove the existence of stable 
checkerboard patterns. 

Unlike the proof in [28], which relies heavily on monotonic- 
ity properties, here we have only assumed that the £2 -gain 
is equal to the dc-gain. 

VI. Conclusions 

In this paper we presented analytical results to predict 
steady-state patterns for large-scale lateral inhibition 
systems. We have shown that equitable partitions provide 
templates for steady state pattern candidates, as they identify 
invariant subspaces where the fate of cells in the same 
class is identical. We proved the existence of steady state 
patterns by relying on the static input-output model of each 
cell and the algebraic properties of the contact graph. One 
limitation in these results is the assumption that the reduced 
graph is bipartite. Therefore, the generalization to a larger 
class of graph partitions, that do not necessarily result in 
bipartite reduced graphs, needs to be investigated. Further 
results will also focus on the case where the cell model 
is multiple-input multiple-output, representing cell-to-cell 
inhibition signals that depend on more than one species. 

Finally, we have analyzed the stability of steady state 
patterns by providing a decomposition into a representative 
subsystem Sb. and a transverse subsystem Sd- We also 
provide a small-gain stability type criterion, which relies only 
on the reduced order subsystem Su to guarantee stability of 
the steady state patterns. 
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